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Biogeochemical  variability  in  the  upper  ocean 

A.  Mahadevan 

Atmospheric  and  Environmental  Research,  Lexington,  MA  - USA 


Abstract.  Differences  in  the  response  times  of  different  biogeochemical 
tracers  to  upper  ocean  forcing  processes  can  lead  to  substantial  differences 
in  their  patterns  of  variability.  An  idealized  tracer  with  a steep  vertical 
concentration  gradient  is  modeled  in  a flow  field  representative  of  an  upper 
ocean  front.  The  surface  distribution  of  the  tracer  is  studied  as  a function  of 
its  response  or  restoration  time  scale. 


Introduction 

Many  quantities,  such  us  temperature,  salinity,  dis- 
solved inorganic  carbon  (DIC),  nitrate,  phosphate,  ra- 
diocarbon and  dissolved  oxygen,  have  strong  vertical 
concentration  gradients  in  the  upper  500  or  so  meters  of 
the  ocean.  These  substances  are  often  forced  in  the  up- 
per ocean,  and  vertical  mixing  and  transport  processes 
that  convey  the  resulting  concentration  changes  to  the 
ocean  interior  are  extremely  small.  Horizontal  advec- 
tion,  on  the  other  hand,  is  relatively  efficient  at  mixing 
and  homogenizing  substances,  and  hence  concentration 
gradients  are  relatively  weak  in  the  horizontal.  Anoma- 
lous concentrations  of  many  biogeochemical  quantities 
in  the  upper  ocean,  are  therefore,  largely  a result  of 
upwelling  that  brings  to  the  surface,  water  with  a sub- 
stantially different  concentration  than  its  surroundings. 
The  substances  may  also  be  altered  by  processes  like 
biological  production,  remineralization,  and  air-sea  ex- 
change. Such  processes  could  induce  horizontal  variabil- 
ity, but  the  variability  induced  by  upwelling  is  typically 
smaller  scaled,  more  intermittent,  and  responsible  for 
most  of  the  variance  in  the  sea-surface  distributions  of 
the  biogeochemical  quantities. 

In  an  earlier  modeling  study  ( Mahadevan  and  Archer, 
2000),  it  was  observed  that  the  nutrient  upwelling  rate 
in  the  pelagic  ocean  is  very  sensitive  to  model  resolu- 
tion. Increasing  the  model  resolution  from  0.4  degrees 
to  0.1  degrees  in  a domain  situated  in  the  interior  of  sub- 
tropical gyre,  resulted  in  a two-  to  three-fold  increase 
in  not  only  the  mean  rate  of  nutrient  supply,  but  also 
in  the  variance  of  the  field.  The  mean  and  variance 
of  other  model  variables  like  temperature  and  pC02  at 
the  sea  surface  are  not  as  sensitive  to  model  resolution 
as  new  production.  The  density  field  in  the  coarse  and 
fine  resolution  model  runs  looks  very  similar  in  its  mean 
and  variance,  but  the  fine  resolution  simulation  shows 
more  folds  in  the  isopycnal  surfaces  and  greater  convo- 
lutedness in  the  frontal  outcrop  formed  by  the  intersec- 


tion of  such  density  surfaces  with  the  euphotic  layer. 
The  transport  of  nutrients  from  the  thermocline  into 
the  mixed  layer  is  primarily  an  isopycnal  process.  Re- 
solving the  length  of  the  front  at  which  the  flow  field 
is  also  highly  strained,  is  thus  extremely  important  for 
capturing  the  vertical  transport  which  occurs  on  scales 
of  the  order  of  10  km  or  less. 

Even  though  all  the  biogeochemical  tracers  in  such  a 
model  are  affected  by  vertical  motion,  the  surface  dis- 
tribution of  some,  for  example  oxygen  and  H2O2,  ap- 
pears more  fine  scaled,  i.e.,  there  is  more  variance  at 
a given  length  scale  than  in  others  like  DIC  and  tem- 
perature. While  there  are  many  reasons  for  differences 
in  the  distributions  of  different  tracers,  I believe  that 
this  is  largely  a result  of  the  difference  in  the  response 
times  of  the  tracers  to  various  forcing  processes.  Sub- 
stances that  respond  quickly  to  upper  ocean  forcing  ap- 
pear in  small  and  intermittent  anomalous  patches  rep- 
resentative of  the  upwelling  scale.  When  substances 
are  slow  to  adjust,  the  anomalous  patches  introduced 
by  upwelling  are  stirred  and  spread  by  horizontal  ad- 
vection.  Hence  the  distribution  of  such  substances  is 
not  as  fine-scaled.  Carbon  dioxide,  for  example,  is  very 
slow  to  equilibrate  with  the  atmosphere  on  account  of 
the  buffering  effect.  While  the  oxygen  in  a 100m  deep 
mixed  layer  may  typically  take  one  month  to  equilibrate 
with  the  atmosphere,  the  DIC  would  take  about  a year. 
The  surface  distribution  of  oxygen  is,  as  a result,  much 
more  patchy  or  intermittent  than  DIC. 

Modeling  Study 

The  change  in  distribution  of  a tracer  c in  the  ocean 
can  be  modeled  using  the  equation 

jt+n -Vc  = 5,  (1) 

where  u is  the  velocity  of  the  fluid  and  5 is  a source 
or  sink  term.  S may  parameterize  processes  like  air- 
sea  exchange  and  biological  production  that  alter  the 
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Figure  1.  Snapshots  of  the  surface  density  field  in  the 
model  shown  on  days  46,  69  and  93  of  the  simulation  help 
visualize  the  evolving  flow  field.  The  flow  is  from  west  to 
east  along  the  the  density  contours. 


quantity  c.  Often  S can  be  expressed  in  terms  of  a 
characteristic  response  time  or  adjustment  time.  In  the 
case  of  salinity  and  temperature,  S could  represent  the 
effect  of  the  heat  or  evaporation-precipitation  flux,  and 
may  be  expressed  as  A_1(c-  Co),  where  c is  the  temper- 
ature or  salinity,  c o is  the  value  to  which  it  is  restored 
at  the  surface,  and  A is  the  restoration  time  scale.  In 
the  case  of  nutrients,  S may  be  written  as  A-1co,  where 
Co  may  be  a combination  of  the  nutrient  and  light,  and 
A a biological  uptake  time.  In  the  case  of  gases  that 
are  subject  to  air-sea  exchange  at  the  ocean  surface,  5 
would  represent  the  gas  flux,  for  which  the  characteris- 
tic time  scale  A is  the  ratio  of  the  mixed  layer  depth  to 
piston  velocity.  Thus  the  quantification  of  the  source  or 
sink  term  S in  terms  of  an  adjustment  time  A is  general 
to  many  tracers  and  processes. 

In  what  follows,  we  use  a three-dimensional  numer- 
ical model  of  an  idealized  ocean  front  to  analyze  the 
surface  distribution  of  a tracer  as  a function  of  its  re- 
sponse time  A.  The  tracer  c is  qualitatively  similar  to 
other  tracers  in  the  ocean;  it  has  a steep  concentra- 
tion gradient  in  the  upper  500  m.  It  is  modeled  using 
equation  (1)  which  is  coupled  to  the  dynamical  model 
described  in  Mahadevan  et  al.  (1996).  The  source  term 
on  the  right  right  hand  side  of  (1)  is  modeled  as 

0,  if  z < zm 

A_1[c-  (l-exp(-^O)],  if  z>zm. 

Here  z is  the  depth  below  the  surface,  zm  can  be  thought 
of  as  a mixed  layer  or  euphotic  layer  depth,  and  zq  is  a 
reference  depth.  In  the  absence  of  advection,  the  tracer 
has  the  steady  state  profile  c = 1 — exp(— (z  — zm)/zo)  at 
depths  greater  than  zm,  and  a value  of  0 at  depths  less 
than  zm.  The  tracer  can  thus  be  thought  of  as  depleted 
in  the  euphotic  layer  (z,  zm),  and  beneath  it,  increasing 
exponentially  with  depth  so  as  to  approach  1 over  an 
e- folding  distance  zq  — zm.  The  euphotic  depth  zm  and 
e-folding  distance  zo  — zm  are  each  taken  to  be  100  m 
in  these  experiments.  Advection  disturbs  this  steady- 
state  distribution  to  which  it  is  constantly  restored  with 
a characteristic  adjustment  time  A.  When  A is  small, 
values  of  the  tracer  anomalous  from  the  steady-state  ex- 
ponential distribution  are  rapidly  annihilated.  When  A 
is  large,  the  anomalous  values  that  are  introduced  pri- 
marily by  vertical  motion,  persist  for  longer,  are  rein- 
forced by  further  upwelling,  and  are  spread  horizontally 
by  horizontal  advection.  Diffusion  is  negligible  at  the 
scales  considered  and  is  not  explicitly  included  in  the 
model. 

The  model  is  set  up  in  a periodic  channel  with  solid 
north-south  boundaries.  The  domain  is  approximately 
256x284  km  in  dimension  and  centered  at  25  °N.  It  is 
initialized  with  lighter  fluid  in  the  upper  100  m of  the 
southern  half  and  denser  fluid  in  the  upper  northern 
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Figure  2.  The  tracer  concentration  in  the  upper  100m 
viewed  on  day  69  of  the  simulation  (middle  panel  of  Fig.  1) 
for  three  different  values  of  the  adjustment  time  A:  2.5,  20 
and  80  days.  The  tracer  distribution  is  more  patchy  and 
the  tracer  concentration  less  for  smaller  values  of  A.  The 
concentration  is  displayed  as  a percentage  of  the  maximum 
concentration,  i.e.  it  is  multiplied  by  100. 


Figure  3.  The  variance  V ( L ) of  the  surface  distribution  of 
tracer  shown  in  Fig.  2 plotted  against  wavelength  (inverse 
of  the  length  scale  L)  on  log-log  axes  for  different  values  of 
A.  V (L)  is  the  average  variance  contained  in  regions  of  area 
L x L.  Smaller  values  of  A result  in  more  variance  within 
a given  length  scale  and  a shallower  curve  in  the  log-log 
space.  The  relationship  between  V{L)  and  A seen  in  this 
distribution  is  qualitatively  the  same  at  other  times  in  the 
simulation. 


half.  Beneath  this  surface  layer,  the  density  in  both 
halves  converges  to  a mean  stratification  that  is  repre- 
sentative of  the  thermocline.  The  model  resolution  is 
0.04°  (approximately  4 km)  in  the  horizontal  and  varies 
in  the  vertical  from  10  m near  the  surface  to  125  m at 
depth.  The  tracer  is  initialized  to  its  steady  state  distri- 
bution and  allowed  to  evolve  with  the  flow  field.  We  use 
five  different  values  of  A:  2.5,  10,  20,  40  and  80  days, 
and  observe  differences  in  the  surface  distributions  of 
the  tracer  as  a function  of  A. 

Results 

The  initial  density  front  gives  rise  to  a geostrophic 
jet  in  the  along-front  direction.  The  jet  and  the  front 
meander  as  the  flow  evolves  giving  rise  to  a more  com- 
plicated eddy  field  that  gradually  decays  (Fig.  1).  The 
surface  distribution  of  the  tracer  that  evolves  with  this 
flow  field  varies  with  A.  It  is  more  patchy  and  has  more 
variance  at  a given  length  scale  for  smaller  values  of  A 
(see  Fig.  2).  The  tracer  distributions  are  analyzed  by 
calculating  the  variance  as  a function  of  length  scale 
[J.W.  Campbell,  personal  communication].  The  vari- 
ance is  computed  over  the  whole  domain,  and  then  over 
subregions  of  size  L,  by  partitioning  and  re-partitioning 
the  original  domain.  V ( L ) is  the  average  variance  over 
all  the  subdomains  of  size  L that  occupy  the  whole  do- 
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Figure  4.  Vertical  profiles  of  the  tracer  at  a late  stage  (day 
93)  of  the  simulation  averaged  over  the  domain  for  different 
values  of  A.  The  upper  ocean  tracer  concentration  increases 
(at  a less  than  linear  rate)  with  increasing  A. 

main.  V(L ) is  normalized  by  the  total  variance  and 
plotted  as  a function  of  the  length  scale  L on  a log-log 
plot  (see  Fig.  3).  We  obtain  shallow  curves  that  have 
small  slope  when  A is  small  and  steeper  curves  for  large 
A.  The  approximate  slope  of  these  curves,  denoted  by 
p,  is  a measure  of  the  intermittency  or  patchiness  of  the 
tracer  distribution.  Small  values  of  p correspond  to  a 
more  intermittent  or  patchy  distribution  with  a large 
amount  of  variance  at  a given  length  scale.  This  kind 
of  distribution  results  from  rapid  adjustment  or  small 
A.  The  converse  is  true  for  large  values  of  p and  A.  The 
interesting  thing  is  that  the  increase  in  p scales  roughly 
as  the  logarithm  of  A. 

The  same  is  true  for  the  average  anomalous  concen- 
tration of  tracer.  Fig.  4 shows  the  mean  vertical  distri- 
bution of  the  tracers  with  different  A towards  the  end  of 
the  simulation.  The  total  amount  of  tracer  in  the  upper 
layer  (z  < zm)  increases  with  A,  but  the  increase  scales 
roughly  as  the  logarithm  of  A for  reasons  that  are  not 
well  understood. 

balance 

Conclusions 

Thus  it  is  shown  that  the  distributions  of  biogeo- 
chemical quantities  in  the  upper  ocean  are  sensitive  to 


their  time  scale  of  response  to  various  forcing  processes. 
This  ought  to  be  taken  into  account  when  trying  to  de- 
sign algorithms  that  quantify  the  surface  distribution  of 
one  biogeochemical  quantity  in  terms  of  another.  The 
patchiness  or  intermittency  of  the  distribution,  quan- 
tified here  in  terms  of  the  approximate  slope  p of  the 
variance  curve  in  log-log  space,  increases  more  or  less 
logarithmically  with  the  adjustment  time  A. 
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